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We calculate the condensate and the vacuum current density induced by external static mag- 
netic fields in (2+l)-dimensions. At the perturbative level, we consider an exponentially decaying 
magnetic field along one cartesian coordinate. Non-perturbatively, we obtain the fermion propa- 
gator in the presence of a uniform magnetic field by solving the Schwinger-Dyson equation in the 
rainbow- ladder approximation. In the large flux limit, we observe that both these quantities, ei- 
ther perturbative (inhomogeneous) and non-perturbative (homogeneous), are proportional to the 
external field, in agreement with early expectations. 
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I. INTRODUCTION 



For massless fermion theories, several aspects of the 
magnetic catalysis of dynamical chiral symmetry break- 
ing, i.e. the formation of a fermion condensate by effects 
of a uniform magnetic field, have been a subject of in- 
tense scrutiny over the past two decades More re- 
cently, the same effect was shown to generate an anoma- 
lous magnetic moment through a dynamically generated 
Bohr magneton inversely proportional to the dynamical 
mass [2j. In this connection, aside from their role as 
toy models in particle physics, theories of fermions in 
(2-|-l)-dimensions have captured the interest of the com- 
munity because the potential applications in condensed 
matter systems, for which the low energy dynamics can 
be described in terms of planar fermions (see, for ex- 
ample 0] and references therein), including graphene in 
the massless version These theories exhibit unique 
features that make them interesting on their own. For 
example, the ground state of an odd number of fermions 
in a uniform magnetic field exhibits a finite value of a 
fermion condensate Q and a parity non-invariant cur- 
rent d, 0] • Parity breaking and gauge noninvariance are 
intimately connected for such a system 0- Gauge in- 
variance can be restored at the expense of introducing a 
parity non-invariant Chern-Simons term in the effective 
action for fermions, or, equivalently, in the action of the 
gauge field. In any case, parity is explicitly broken. The 
formation of the condensate by effects of a magnetic field 
is important, for example, for planar doped antiferromag- 



* Electronic address: raya@ ifm.umich.mx| 

t Electronic address: edward . reyes @n ucleares . unam. mx| 
■tAlso at: Institute de Fisica y Matematicas, Universidad Michoa- 
cana de San Nicolas de Hidalgo. 



nets. Such systems are relevant to the physics of high- 
Tc superconductors [1, |1] in the surface region of these 
materials, where, due to the Meissner effect, the exter- 
nal magnetic field can penetrate the sample. In Ref. Q, 
the dynamical mass gap generated by an intense external 
uniform magnetic field was studied in the reducible for- 
mulation of parity-invariant quantum electrodynamics in 
the plane, QED3, solving the Schwinger-Dyson equation 
in the rainbow-ladder truncation and constant-mass ap- 
proximation, in a dimensionally-reduced variant of the 
well know expression of magnetic catalysis in (3+1)- 
dimensions. In that work, the proper-time representa- 
tion [lo| of the fermion propagator was used. The dy- 
namical formation of this gap was found to be connected 
to the enhancement of the superconducting gap in the 
strong- [/ Hubbard model of Ref. Q. It is important 
to point out that the magnetic field in the surface re- 
gion of these materials is not uniformly distributed in 
space. Therefore, it becomes important to consider spa- 
tial anisotropics in the formation of condensates. In this 
connection, through the second quantized solutions to 
the parity-invariant Dirac equation in a reducible rep- 
resentation of the 7^-matrices, it was shown [ll| that 
for non-uniform static magnetic fields of intense flux, the 
inhomogeneous condensate follows the spatial profile of 
the applied field, in a local version of the Aharonov- 
Casher integrated relation [l3|, which is essentially the 
Landau's degeneracy-flux relation [l^. Such a relation is 
connected with mathematical index theorems [11]. The 
inhomogeneous flelds considered in Ref. [llj include fields 
which vary along the radial or a cartesian coordinate. 
The contribution to the condensate in the massless limit 
was found to come only from the threshold states, for 
which the mass equals in magnitude its energy. 

In this work, we are interested in obtaining the fermion 
condensate and vacuum current density induced by ho- 
mogeneous and inhomogeneous magnetic fields, but our 
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approach is different. We derive tfiese quantities directly 
from the fermion propagator. For this purpose, we use 
the so-cahed Ritus eigenfunctions method [T^, which 
provides an elegant and powerful method to diagonalize 
the fermion Green's function in momentum space in the 
presence of external fields. The method was originally 
developed for spin- 1/2 fermions, and was later extended 
to the case of charged bosons [16i]. In both these cases, 
uniform fields were considered. Here, we present a gen- 
eralization of the Ritus approach to the case of inhomo- 
geneous fields in (2-|-l)-dimensions [13 • We assume the 
field aligned perpendicularly to the plane of motion of 
fermions. The diagonalization of the Green's function in 
momentum space for uniform and non-uniform magnetic 
fields might also be useful for the theoretical studies de- 
rived from the measurement of half integer quantum Hall 
effect and zero-energy Landau level in graphene at room 
temperature [l^. Moreover, the method allows to write 
down general expressions for the condensate and the in- 
duced current density presented in Eqs. (|1T|) - (H51) . We 
consider two special field configurations. First we im- 
plement Ritus method to express the tree-level fermion 
propagator in the presence of an exponentially damped 
magnetic field along one spatial dimension. Such a field 
occurs, to a good approximation, inside the London pen- 
etration depth of a type-I superconductor if a semicon- 
ductor heterostructure -with narrow quantum well- is in- 
troduced perpendicularly to the planar surface of the 
superconductor into a narrow slit, and a homogeneous 
magnetic field is applied parallel to the surface of the su- 
perconductor [19|. As a second example, we obtain the 
non-perturbative fermion propagator in QED3 with the 
same assumptions that Ref. [9] , but expanding the prop- 
agator in the Ritus eigenfunctions. From this propagator 
we derive the non-perturbative condensate and induced 
charge density. This article is organized as follows: in 
Sect.ini we describe the Lagrangian for Dirac fermions in 
(2-1-1 )-dimensions and in the presence of external fields, 
with emphasis in the symmetry properties of the mass 
terms. Section IlIII is devoted to the fermion propagator 
in the presence of an inliomogeneous magnetic field which 
is exponentially damped along one spatial direction [l3| 
within the Ritus formalism [l5|. The perturbative con- 
densate and induced vacuum current density for this field 
are discussed. The non-perturbative propagator in a uni- 
form magnetic field and the corresponding homogeneous 
condensate and vacuum charge density are discussed in 
Sect. IIVI Final remarks are presented in Sect. |V1 

II. PLANAR FERMIONS 

We start from the Dirac Lagrangian in external elec- 
tromagnetic fields 

£0 = 1^(7 • n - m)V', (1) 

where H^ = id^ + eA^. For a more detailed presen- 
tation of the symmetries of this Lagrangian, see, for 



instance, Ref. [1^. In (2-|-l)-dimensions, only three 
Dirac matrices are required to fulfill the Clifford algebra 
{7^77''} — '^g^" ■ The lowest dimensional representation 
of these matrices is 2 x 2 and hence we can choose them 
to be proportional to the Pauli matrices as 

7° = f^3, 7^=«o-ii 7^ = icr2- (2) 

We call this representation A. In this representation, 
fermions posses only one spin orientation. Additionally, 
there exists a second inequivalent representation, labeled 
B, which can be chosen as follows 

7° = era, 7^=^0-1, 72 = -10-2, (3) 

in which fermions have the opposite spin orientation 
alone. In graphene, representations A and B are required 
to describe two different species of massless fermions in 
each triangular sub-lattice of the honeycomb lattice [2l[ . 
Chiral symmetry cannot be defined for either A or B, 
because there is no 2 x 2 matrix analogous to '75'. More- 
over, the mass term rmpip in the Lagrangian is non- 
invariant under the parity transformation xi — > —xi and 
Ai{xi,X2) —Ai{—xi,X2) for these representations. 

The two species with their respective spin orientations 
can be conveniently combined into a 4-component spinor 
with a 4 X 4 representation of the Dirac matrices, say, 

for k — 1,2, which we label C, and 

7^=^(;;), 7^=nV7v^. (_;;). (5) 

Here I is the identity matrix (we use the same symbol 
in any dimensionality). In this reducible representation, 
miptp is parity invariant. Furthermore, in the massless 
limit, the Lagrangian ([1]) possesses a global U{2) fla- 
vor symmetry, with generators I, 7"^, 7^ and [7'^, 7^], 
corresponding to the interchange of the two irreducible 
species. The ordinary mass term breaks this symmetry. 
However, there exists a second mass term of the form 
nio'tpTilj, with T = [7'',7^]/2 — diag(I, —I), which in con- 
densed matter literature is often referred to as Haldane 
mass term [2^ . This mass term is invariant under fla- 
vor symmetry, but breaks parity. In this case, the parity 
non-invariant Dirac Lagrangian takes the form 

= "0(7 • H — TO — moT)ip. (6) 

In order to explicitly separate the physical fermion con- 
tent of this Lagrangian, we introduce the chiral-like pro- 
jectors x± = (Ii''")/2 and define the "right-handed" ^+ 
and "left-handed" V- fields as '0± = x±V'- Then, the 
Lagrangian acquires the form 

£d = '0+(7 • n - m+)xp+ -I- '0_(7 • H - m^)ip^, (7) 

where m± = m±mo. Thus, in this form, the Lagrangian 
is neatly seen to describe two different fermion species. 
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and the effect of the parity-violating mass term is to re- 
move the mass degeneracy between them. Below we ob- 
tain the fermion propagator for Lagrangians ([T]) and ([7]) 
in external magnetic fields. 



III. PROPAGATOR IN INHOMOGENEOUS 
MAGNETIC FIELDS 

In this section we obtain the fermion propagator in an 
inhomogeneous magnetic field perpendicular to the plane 
of motion of the fermions in (2-|-l)-dimensions within the 
Ritus formalism [15| . Working in a Landau-like gauge, we 
choose Af^ = (0, 0, W{x)) such that the profile of the field, 
which we consider varying along the x-axis, is B(x) — 
W'{x) = dxW{x). Details of this derivation can be found 
in Ref. [13] ■ The Green's function for the Dirac equation, 
S{z, z'), satisfies 



(7-n-m)S'(z,z') =<5(3)(z-z'). 



(8) 



with z = {t,x,y). Since S{z,z') commutes with (7 -H)^, 
we expand it on the basis of its eigenfunctions. 



(7-n)X(z)=pX(z). 



(9) 



Without loss of generality, we work only with the ir- 
reducible representation A, and specify how results are 
modified in other representations. In this case, (7-11)^ = 
lP+ea3W'{x), and thus the Ritus eigenfunctions become 

(10) 

With the aid of the property 7^7'^ — gf^"^ — ie^'^^'jx, these 
can be expressed in the more convenient form 

E(^\z) - i?p,+i(z)A(+) + i?p,_i(z)A(-) , (11) 

with 

I ± ^7^72 



A(±) 



(12) 



being the spin projectors "5]. In the above expressions, 
the subscript p — {po,P2, k) denotes the eigenvalues of the 
operators idt, —idy, and H ~ — (7-n)^+nQ, respectively, 
a — ±1 are the eigenvalues of 173 and 

E,.A^) = N,e-^^P^''P^y^F^,^^{x), (13) 

with No- a normalization constant. F^p^{x) satisfies 

[dl - {-P2 + eWix))^ + eaW'ix) + k] F^^^ (x) = 0, 

(14) 

which is the equation of the Pauli Hamiltonian with con- 
strained vector potential, mass to = 1/2 and gyromag- 
netic factor g — 2. This Hamiltonian turns out to be 
supersymmetric in the sense of SUSY-quantum mechan- 
ics [23]. F^p^{x) and F^^^{x) are the solutions of the 
respective StJSY-partner potentials 



Regrettably, the solution to Eq. (fT4)) for arbitrary W{x) 
is unknown. For certain inhomogeneous fields which are 
translationally invariant along one spatial direction, the 
solutions can be expressed in terms of an orthogonal sys- 
tem of functions [23|. However, the following important 
property of the Ritus eigenfunctions. 



(7-H)Ep(z)-Ep(z)(7.p), 



(16) 



where the three-momentum vector p satisfies p^ — p^ — 
Pq — k with Pfj^ — {pQ,0,\/k) is valid in the general 
case [17| . The explicit matrix form of Eq. ([T6| is 



idtEp^+i{z) D-Ep-i{z) 
D+Ep^+i{z) -idtEp^^i{z) 

^ ( PoEp^+i{z) 
\ VkEp^^iiz) 



-VkEp^+i{z) 
-PoEp,^i{z) 



,(17) 



with D± — —dx ± {—'idy — eW{x)). The relations follow- 
ing from diagonal components of this matrix equation 
are directly inferred from the properties of the Ep^a{z) 
functions, while the off-diagonal components can be cast 
in the form of the 'kinetic-balance' system of equations 



D^Ep^^iiz) = -VkEp^+i{z), 
D+Ep^+i{z) = VkEp^.i{z). 



(18) 



Thus the Ep_„{z) functions which satisfy these expres- 
sions also satisfy the identity (|16p . Its explicit form, 
however, will depend on the field under consideration. 



A. Exponentially Decaying Magnetic Field 

In order to proceed further, we concentrate the discus- 
sion to the case of an exponentially decaying magnetic 
field B{x) = i?e~"^. Solutions to the Dirac equation in 
this field has been studied in Ref. [l^ [13, [l^- Such 
a field occurs inside the penetration depth of a type- 
I superconductor when a semiconductor heterostructure 
with narrow quantum well is introduced perpendicularly 
to the planar surface of the superconductor, and a ho- 
mogeneous magnetic field is applied parallel to this sur- 
face [l^. An exponentially decaying field can be de- 
scribed by W{x) = — (i?/d)[exp{— dx}-l]. This choice 
allows to directly recover the uniform magnetic field case 
of Sect. Iiy I b v setting a = 0. Defining the dimensionless 
variables Il9l IHi 



eB 



1 



[ai{x) 



2 ' 



and 



a 0? [ci£(a;o)] 



(19) 



(20) 



V±(a;) = (-P2 + eM{x)f ± eW'{x). 



(15) 



where l{x) is the magnetic length, ^(xo) its local counter- 
part [l^,[2j andp2 = pi^&B ja, the Pauli equation 
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takes the form 
Id 



(2s 



(21) 

Here, = kjo?. The behavior at small and large ^ 
of this equation suggests Fj^p^iO ~ f ^e~^a;'^(^), with 
= — Inserting this ansatz, Eq. ((2T|) becomes 



c^^(^.) = 



(22) 



with A = 2/3 + 1 and 77^ = -s + /3 + 1/2 - cr/2 
and g — 2^. The general solution to this equation is 
= C(i) iFi[77^, A, £i] +C(2)U[T]^,\,g], where iFi[a,6,a;] 
and C/[a, 6, x] are, respectively, confluent hypergeometric 
functions of the first and second kind. The irregular be- 
havior of U[a,b,x — 0] enforces C(2) = 0. Moreover, to 
preserve the asymptotic behavior, it is required that rj^r 
is a negative integer, — |77o.|=0,l,2,... This implies that 



(23) 



The energy eigenvalues for a particle on-shell, = m^, 
are conveniently written as 



kn=p\- {P2 - nay 



(24) 



Notice that fc„ explicitly depends on the momentum p2 , 
unlike the case of the uniform magnetic field case, even in 
(3+l)-dimensions. For such a particle, the energy Pq = 
kn + m? cannot be tachyonic. This fact restricts (3 = 
s — n > 0, which can be achieved so long as p2 > an. 
We observe that for fc+i, we can write — jry+ij = n — 
1, and for /c-i, — |f7-i| ~ n. Regarding the solutions 
uj'^{g), for n = 0, ujQ^{g) does not exist, whereas ujQ^{g) 
is nondegenerate. With all the above, defining ui!^_i{g) — 
0, the solutions to Eq. acquire the form 



From the identity 



l,g 



(25) 



Ll,{x)=^-^^^ ,F,[-m,K+\,x], (26) 

where LJ^(a;) are the associated Laguerre polynomials 
and (k)„ — T{k + m)/T{K) is the Pochhammer symbol, 
solutions to Eq. (^11) can be expressed as 25] 



where 



"22 



(27) 



(28) 



are the Laguerre functions |26l | , which form a complete set 
of orthogonal functions in the interval (0,cx)) and verify 
the orthogonality relation 

dx r^{x)r^, (x) = nm+J + l) ^29) 
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FIG. 1: SUSY-parter potentials H15|) and energy eigenvalues 
for the exponentially decaying magnetic field for s = 5. The 
number of bound states grows linearly with s. The scale of 
the plot is set by a = eB — 1. 



This relation is important to understand the spectrum of 
bound states of Eq. (PT|) below. 

Now, the energy eigenvalues in Eq. ([M)) suggest the 
identification of the quantum number n with the Landau 
level index. This is indeed the case, because such eigen- 
values, in the limit d 0, reduce to the Landau levels 
for a uniform magnetic field. 



wl{xo)'2n , 



(30) 



where Wc{xo) — l/£{xo) is the local cyclotron frequency. 
Moreover, because s is a real parameter, the number of 
normalizable solutions for a given s, restricts the Landau 
level index n to take the values n = 0, 1,2, ...,<;— 1 where 
— [s] is the integer part of s. This is evident from Fig.[Tl 
where the SUSY-partner potentials (ITSt are shown as a 
function of x for fixed s = 5. The energy eigenvalues fc„ 
for the normalizable states are also displayed. 



With all the above, the solutions Ep^cr{z) in Eq. ([13 
acquire the explicit forms 



Ep,+i{z) 



a f 2n!(s — n) 



1/2 



-ipot+ip2y 



2tt VF(2s - n + 1)^ 

xe-^/^g'^'-^^Ll^'-'^HQ), 



2-K \ F(2s - n) 



(31) 



Inserting these expressions into Eq. ((TU)) . we obtain the 
Ritus eigenfunctions for an exponentially decaying static 
magnetic field. It is straightforward to check that these 
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P 



verify the orthogonality and closure relations 

Ep{z)Ep{z') = (5(3) (z - z'), (32) 
where Ep = j"EIj° , 

6'^^\p-p') = Sn,n'5s,s'S{po - p'o), (33) 

and the projector 0, Hzl 

n(n) A(+)5„,o+II(l-5n,o), (34) 

determines only one spin projection in the lowest Lan- 
dau level (LLL). The symbol J^d^p indicates that the 
integration might represent a sum, depending upon the 
continuous or discrete nature of the components oi p = 
{po,P2,k). In our example, 



(35) 



Equations ([5^ follow directly from the complete and or- 
thogonal character of the solutions 



B. Condensate and Induced Electric Current 

Physically, Ritus eigenfunctions Ep{z) correspond to 
the asymptotic states of electrons with momentum p in 
the external field. Therefore, we can use these functions 
to expand S{z^ z') in momentum space in the same way 
plane waves are used to define the Fourier transform. 



d-'p d'p' Ep{z)S{p,p')Ep,{z') . (36) 



Inserting this Green's function in Eq. (jS]), using the prop- 
erties (jl6p and (j32p . the propagator in momentum space 
takes the form 



where 



S{p,p') = PHp-p')Ilin)Sip), 



S{p)^ 



^ ■ p — m 



(37) 



(38) 



with p — (po,0, \/p2 ~ (P2 — naY). This simple form 
allows a direct calculation of physical observables, like 
the ferniion condensate 

{^^)^Tr{iS{z,z)}, (39) 

and the induced vacuum current density 

= -ieTr{7^S'(z,z)}. (40) 



These acquire the general form 



X [|£;p,+i(z)|2 + |i?p,_i(z)|2], (41) 
-ie^^d^p ™ 

X [\EpM^)\^-\E^..-i{^)?], (42) 

0, fc = l,2. (43) 



J A 



J A 



We want to emphasize that these expressions are valid 
for any profile of the magnetic field so long as we know 
the solutions to the Pauli equation (ITil) . from which we 
can build the functions Ep^„{z) of Eq. p^ . 

Inserting the explicit solutions, the condensate and 
charge density are 



27r 



Y — 



k?=0 

oo oo 



|m| Vr(2? + i) 

oo 

E E 



„=i<;=„+i V"^(2<rn-n2)-f m2 

r(2<7- (n- i))y 

r(2,-n) ) ^^(-1) ^^^^ 



(44) 



2tt 



OO ^ 

Y — 



\m\ Vr(2^ + i) 



OO oo 



( n!(<;-n) \ , ,2i,-n) ^^^^ 



\T{2,-{n-l)) '^^^^ 

(n-l)!(^-n) \ ,j2{,-n),..2 
r(2,-n) ) ^^(-1) ^^^J 



(45) 



Here we have only integrated over nornializable states 
and separated explicitly the contribution from the LLL. 
Now, because ^ takes only discrete values, then Landau 
levels are highly degenerated, except the LLL. Therefore, 
for fields of moderated flux, we need to perform and regu- 
larize the remaining sums. However, in the case of intense 
field, the leading contribution comes from the LLL. So, 
in this regime, as m — )■ 0, 



3 A = ^sinh(£()sgn(TO), (46) 



47r 



ge ^ sinh(p)sgn(m) 



(47) 
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In the intense flux limit, the above expressions become 

o2 



in 



Be~°"'sgn{m), (48) 



mT°'^=' = -Be-"-sgn(m), (49) 

and comprise the main results of this section. 

Translation of these findings to representation B of 
Dirac matrices given in Eq. ([3]) is straightforward. Ri- 
tus eigenfunctions can be constructed as 



(50) 



and then, we find that {ipip)B — ('0^)yt: 3b ~ ^^.d 
A — Jb — ^ fo'^ ^ — whereas for the reducible 
representation C, Eq. the Ritus eigenfunctions are 



Ef (z) 



i-^\z) 







(51) 



The block-diagonal structure emphasizes the existence of 
two fermion species, each with a different mass. Thus, 
formally, (■0V)c = ('0'0)>t(™+) + (V''0)>t(»7i-) and = 
j^(TO_i_)-l-j'^(r7i_). At the end, we take the massless limit 
m± — >■ of these expressions. 



IV. PROPAGATOR IN UNIFORM MAGNETIC 
FIELDS 

In this section, we obtain the fermion propagator in 
a uniform magnetic field by solving the corresponding 
Schwinger-Dyson equation (SDE) in QED3. As a first 
step, we construct the Ritus eigenfunctions with the pro- 
cedure outlined earlier. A uniform magnetic field can 
be specified by the choice W{x) = Bx. Then, we sim- 
plify the Pauli equation ((T4)) replacing k — >■ 2\eB\k and 
making the change of variable -q — \j2\eB\[x — p2/{eB)], 
obtaining 



,21 



^ + ^-+^sgn(eB)-^ 



KpM^O- (52) 



Solutions are parabolic cylinder functions Z3„ {x) of order 
n = k + (TSgn(ei3)/2 — 1/2. The normalized Ep^a- are 



(7r|eB|)4 
27r3/2fc!^ 
(7r|ei?|) 



27r3/2(fc- 1)! 



From these, we can build up the Ritus eigenfunctions 
in A representation, Eq. (ITU|) . We use these functions 
to solve the Schwinger-Dyson equation for the fermion 
propagator in the rainbow-ladder approximation Q. 



To this end, let us recall that the full fermion propa- 
gator verifies 

[7 • n - z')]G{z, z') = (5(3) (z - z'), (54) 

where Y,{z,z') is the fermion self-energy. In Ritus for- 
malism, the full propagator is expressed as 

G{z,z') ^^d^P d^p' Ep{z)G{p,p')Ep>{z') , (55) 



with 



and 



where 



Gip,p')^P\p-p')Il{n)G{p), 



G{p) 



1 



7-p-S(j3)' 



t{p)^^-p Z{p) + M{p). 



(56) 



(57) 



(58) 



On the other hand, in the rainbow-ladder approximation, 
the self-energy takes the form 



I](z, z') = -^6^7^6(2, z')-i''D^,{z - z'), (59) 



where 



D^^(z - z') 



(27r)3 g2 _ i£ 



x(ff,. + (A-i)M;^ 



(60) 



is the bare photon propagator and A is the covariant 
gauge fixing parameter. Ritus eigenfunctions allow to 
diagonalize the self-energy in momentum space as 



S(p,p') = J d^zd^z' Ep{z)^{z,z')Ep'{z') 

= S^^^ip~p')U{n)±{p). (61) 
Combining these ingredients, Eq. ()59p becomes 

d^zd^z' Ep(z)S(z,z')Ep'(z') = 

-ie^ / d^zd^z'D^^iz - z')Ep{z)-/'' x 



P 



d^p d^p' Ep{x)G{p,p')Ep,(z') 



j-Ep,{z') .(62) 



In Feynman gauge, Z{p) — 1. Moreover, we consider 
the LLL approximation, and neglect the momentum de- 
pendence of the self-energy in the SDE equation, i.e., 
we assume the so-called constant mass approximation, 
M(0) = rriiiyn. In doing so, we only take into account 
the part of the self-energy which is proportional to the 
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FIG. 2: Dynamical mass in the constant mass approximation, 
Eq. (|66p as a function of a and eB. 



identity matrix. There is no general principle which guar- 
antees the validity of this approximation. However, in 
in Ref. [9|, this approximation has been established to 
be relieble in QED3 through the Schwinger proper-time 
approach. Here we work under the same assumptions 
within Ritus formalism. After a lengthy but standard 
procedure, the SDE equation reduces to 



1 = le^^neB 



q e' 



1 



(27r)3 



(63) 



''dyn 



where = Qq + q^, = ql + ?! and Q = Q/ \j2eB for 
Q = qo, qi, q2. Defining a ~ e^/(47r), after straightfor- 
ward integration, the above expression reduces to 



1 = 



^dyn 



zTrerf 



Ei 



'"-dyn 

2eB 



(64) 

where —ieviiix) is the error function with complex argu- 
ment and Ei(a;) is the exponential integral function. For 
consistency of the approximation, we require eB ^ m, in 
such a way that mdym obeys the transcendental relation 



1 = 



a 



'^dyji 



■log 



2eBe-^'' 



dyn 



(65) 



with ~ 0.577216 being the Euler constant. Thus, 



nidyn = '2.aW 



_ TB 

e 2 



V2Ib\ 



(66) 



where W(x) is the Lambert W function, i.e., the inverse 
of the function j(w) = we^ for any complex number w. 
In Fig. [5] we display rridyn as a function of a and eB. It 
is positive definite. The result in Eq. (j66p was derived in 
Rcf. 9,] 



Inserting the non-perturbative propagator 
Eq. wc obtain the charge density 



.Q e uidynB 



dpo 



Pi 



dyn 



e^B 
' 47r ' 



into 



(67) 



in agreement with well known perturbative results Q 
identifying m = nidyn- Notice that in Eq. (j42|) . the dif- 
ference [|£:p, + i(z)|2- |^p,_i(z)|2] is such that there ex- 
ists a neat cancellation of the contribution to the charge 
density between subsequent Landau levels, and only the 
LLL contribution prevails. 

The fermion condensate, on the other hand, becomes 



n=o,<j=i _ 

" 47r 



(68) 



Equivalent expressions for this quantity have been ob- 
tained by different methods [5J after the same identifica- 
tion m — nidyn- 



V. FINAL REMARKS 

In this article we have studied the formation of conden- 
sates and vacuum electric current densities of the ground 
state of massless fermions in (2-|-l)-dimensions by ho- 
mogeneous and inhomogeneous magnetic fields. These 
quantities were extracted directly from the fermion prop- 
agator. The effects of an external magnetic field were 
included within the Ritus eigenfunctions approach fis} . 
which was generalized to incorporate magnetic fields of 
arbitrary spatial profile. The class of field configurations 
that can be considered within this formalism are those 
for which Eq. ([T4l) can be solved p/Tj], and similar con- 
clusions are expected in all these cases General ex- 
pressions for (V'V') and are presented in Eqs. (HT|) - (H5| . 
Although we have worked out explicitly the derivation of 
these quantities only in the irreducible representation of 
the Dirac matrices, Eq. ensuring that only one spin 
orientation for fermions enter in the LLL, we have spec- 
ified how our findings can be translated to the second 
inequivalent ([3]) and the reducible dH representations, 
where we have also considered parity non-invariant mass 
terms. 

In the large magnetic flux regime, we have seen that 
both (ip'ip) and are proportional to the external field. 
In perturbation theory, the local relation between the in- 
homogeneous condensate and the fiux can be interpreted 
as a local form of the Aharonov-Casher relation 12|, as 
anticipated earlier in The induced current density, 
in turn, has the form 



.fix) 



(69) 



with *F^'{x) = e^"'^F^x{x)/2. Therefore, it is gauge 
invariant and conserved. We observe that there ex- 
ists a LLL dominance for the formation of the conden- 
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sate and charge density for intense inhomogeneous mag- 
netic fields. For the non-perturbative uniform conden- 
sate and induced current results, an interesting question 
that naturally arises is whether in the magnetic catal- 
ysis scenario in (2-|-l)-dimensions a Chern-Simons term 
of non-perturbative origin should be considered to cure 
the anomaly. However, for uniform fields, such a term 
vanishes formally for our choice of the vector potential 
^/x 0- Thus the computation of the induced current 
alone is not sufficient to deduce the presence of such a 
term in the complete effective action. The effective action 
r for the gauge field has to be computed simultaneously 
and then, the Chern-Simons term can be inferred from [7] 



This and other additional effects of homogeneous and 
inhomogeneous magnetic fields, like the dynamical gen- 
eration of mass and anomalous magnetic moment , for 
the Lagrangian ([7]) are currently being considered and 
will be presented elsewhere. 
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